from jsmodel_bistable_switch_server import *
from pylab import plot
from pylab import show
from pylab import legend
from pylab import figure
from pylab import semilogy


t = 9500
p1 = 191
p2 = 190
p3 = 190
p4 = 190
p5 = 200
p6 = 200
pcs = [3500]
pchl = 0
pcd = [500]
rgs = [0]
rghl = 0
rgd = [500]
thresh = 1e2
req = False
#Hill coefficient: inducer 1, inducer2, tal97, tal57
hk = [1.0,1.0,3.0,3.0]
switch = BistableSwitch(True, t, p1, p2, p3, p4, p5, p6, req, pcs,  pchl, pcd, rgs, rghl, rgd, hk, thresh)
figs = switch.run()
f0 = figs[0]
f1 = figs[1]
f2 = figs[2]
f3 = figs[3]
plot()
show()
